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Abstract 

We find that in two dimensions wires can crumple into different morphologies and present the 
associated morphological phase diagram. Our results are based on experiments with different 
metallic wires and confirmed by numerical simulations using a discrete element model. We show 
that during crumpling, the number of loops increases according to a power-law with different 
exponents in each morphology. Furthermore, we observe a power-law divergence of the structure's 
bulk stiffness similar to what is observed in forced crumpling of membranes. 

PACS numbers: 05.45.-a, 68.35.Rh, 89.75.Da 
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Crumpling is omnipresent in nature, occurring on all length scales, ranging from blood 
cells that crumple in order to pass through capillaries up to the formation of the Swiss 
Alps. The crumpling of spatially extended membranes became a challenging research topic 
of strong interest over the last decades. Thorough research, both experimental, theoretical 



and numerical was conducted to discover basic statistical properties 



of strength and energy [3|, |4|, |5| and landscapes in phase space 
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In this letter, we present the first morphological phase diagram of crumpled wires in 
two dimensions. Metal wires with different material behavior are crumpled by injecting 
them into a circular cavity from opposing sides. To study the material dependence of the 
morphology, we keep the general geometric setup fixed, instead of allowing e.g. for the 
injection angle or cavity shape to be varied like in 8j. Our experiments show that plastic 
yielding and friction are the essential parameters determining the morphology. We find 
excellent agreement with computer simulations and construct a full morphological phase 
diagram. We further analyze the morphologies by means of the number of loops as function 
of the packing density, which exhibits a power-law behavior with different exponents for 
each morphology. Finally we show that the stiffness of the crumpled structures follows a 
power law increase in the packing fraction, similar to observations on membranes 4, 
The stiffness, however, does not fully diverge due to the finite wire thickness, resulting in 
jammed states where percolating force chains are formed by contacting wire segments. 

In our experiments a wire is pushed from opposing sides into a metallic cylindrical cavity 
of radius r=10cm having a width of one wire diameter d with a transparent acrylic glass top 
to allow for observation. The setup is based on experiments by Gomes [§], with the crucial 
difference, that we are pushing the wire in a controlled way by two sets of counter rotating 
rollers driven by large forces. With this setup we crumple wires of steel (d=0.8mm) and 
brazing solder (d=lmm) - materials with differing mechanical properties. 

We model the wire by point masses, connected by tensile springs. Bending stiffness is 
considered by rotational springs attached to each node. We choose Hook's law 

Fi = —ki(X - Aj) and M { = -k 2 Q (1) 
with the rest length A , the angle Q at node i and the two proportionality constants k\ 
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FIG. 1: Crumpled wires in a circular cavity. Experiments (top row) and simulations (bottom row), 
showing the three different morphologies. Materials are steel wire (1.5125) d=0.8mm without 
and with silicon oil (left) and (middle) and brazing solder (Sn99Cul) d=lmm (right). Simulation 
parameters are /=100, K0=O.O4;^ s j=5;s=O.O8 for the classic (left), Ke=7r/2;/i s j=0;s=0.04 for the 
spiral (middle) and Ke=0.02;/i st =5;s=0.04 for the plastic (right) morphology. Simulation snapshots 
are shown for similar packing densities as in the respective experiments. 



and k 2 for the tensile forces and bending moments respectively. From the continuum limit 



ll| it follows that k± = EA/X and k 2 = EI/X with Young's modulus E, wire cross- 



section A and second moment of inertia I. For the numerical results we fixed £7=1275, 
corresponding to /ci=1000 and k 2 =60 for a circular wire with diameter d—1 and Ao=l. We 
use the dimensionless ratio / = r/d of cavity radius to wire diameter to specify the effective 
system size. 

We describe plastic deformation of wires in the rotational springs, using a simple linear 
flow rule with yield threshold moment of k 2 ag an d slope sk 2 (0<s <1). Therefore Mj in 
EqJT] holds for ( < ag, while for £ > ag M = —k 2 a e — (£ — ag)sk 2 . Unloading is always 
along the elastic path with slope k 2 . This is a linear approximation to the bending stress- 
strain relation, as can be found for example in Refs. [12]. In the following, we will use the 
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discretization-independent yield curvature Kg — a$/ Ao as plasticity parameter. 

For a realistic simulation considering friction at wire-wire and wall-wire contacts proved 
to be crucial. We therefore implemented a simple Coulomb's friction law F st < /i s t-F/v with 
friction coefficient fi st and force normal to the contacts F^. Above fist^N, dynamic friction 
sets in with = /i^-F/v, opposed to the relative tangential movement of contacts. Below 
a threshold velocity Vth, static friction sets in again 13J. For numerical reasons, we add a 
small viscous damping on the translational and rotational degrees of freedom. 

We follow the time evolution of the system by integrating the equations of motion of all 
nodes. We choose /=100 to match the experimental setup. During the simulation we push 
new elements into the cavity, starting from an initially cavity-spanning straight wire. The 
simulation stops, when a sudden increase of the injection force appears. At this point, we 
consider the system as jammed. The experiments end, when the resistance of the crumpled 
structure becomes so large, that the rollers are not able to push more wire into the cavity. 

The crumpling starts with an initial buckling of the cavity-spanning wire in the up- or 
downwards direction. The symmetry of the wire is broken by small deviations from the 
ideal positions. This solution is stable until the wire contacts the cavity wall for the first 
time. Depending on yield threshold Kg and static wire friction coefficient fi st , three different 
morphologies can be distinguished: 

(1) With high cavity friction fi st a morphology emerges that we call classic, since it 
corresponds to the observation of Ref. 8fl. As can be seen in Fig{T](left), the wire immediately 
forms cascades of loops of decreasing sizes. Since wire is inserted on both sides, strong 
symmetry is present. In simulations of purely elastic wires, this symmetry is preserved until 
the structure is jammed, while experimentally, it is broken when plastic flow sets in for 
high packing densities. Cascading loops and strong symmetry are the two criteria for this 
classical morphology. 

(2) Elastic wires with low fi st form a spiral pattern, with the turning direction chosen 
spontaneously (see FigJITmiddle)) comparable to Ref.[9(. The spiral winds up more when 
more wire is inserted. Simulations reveal that this pattern is maintained for purely elastic 
wires until the entire cavity is filled. In the experiment, however, a critical packing density 
might exist where plastic deformations set in. In that case, loop cascades appear on the left 
and right side of the spiral, i.e. a mixing of the classic and the spiral morphology can be 
observed. 
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(3) Typical for the plastic morphology is the loss of symmetry. This third morphology 
arises for highly plastic wires (FigJT](right)). It starts similarly as the classic one by form- 
ing cascades of loops. In contrast to the classical morphology the axial symmetry is soon 
broken and large rearrangements of existing structures occur. Such reordering is observed 
on length scales ranging from a single loop up to the rearrangements of entire cascades. As 
a consequence, higher packing densities can be obtained than for the other morphologies. 
Note that the plastic phase is defined by the disorder and not by the material. Plastic and 
classic morphologies are separated by comparing the spatial distribution of curvature, which 
is concentrated near the cavity border for the classic phase (cf. the location of loops in Fig. 




Experimentally, we produced the classic, spiral and plastic morphologies using steel wire 
without and with silicon oil or brazing solder, each one represented by one point in the phase 
space. We investigated this phase space numerically in detail by changing yield curvature Kg 
from to 0.06 and friction fi st from to 0.6, resulting in the morphological phase diagram 
of FigfSJ Simulations are limited to /=50, to be able to present a precise phase diagram 
in reasonable CPU-time, however preliminary results for /=100 do not show significant 
differences. For Kg = we allow plastic deformation right from the start, while Kg = 7r/A 
represents the elastic case. For /i st >0.6 and Kg>0.06 we observed no further influence on the 
morphology. For all simulations fci=1000, /c2=60 and s=0.05 were kept constant. We find 
that a direct transition from the spiral to the classic phase is not possible and a reentrant 
phenomenon is observed at 0.2</i st <0.35 with increasing Kg. For an explanation, consider 
the two requirements for building the symmetric structures of the classic phase: First, 
friction stabilizes the structures during packing. With pL st between 0.08 and 0.18, this 
constraint is not sufficient and the loops rotate and rearrange. Consequently, the plastic 
phase emerges between the spiral and classic one. However, also plastic deformation has a 
stabilizing effect by dissipating stored energy, limiting the system's capability to rearrange. 
To find the classic phase at 0.2<yU si <0.35, small Kg is required to prevent rearrangements, 
along with friction. For too small Kg, no cavity spanning loops form and the classic phase 
is unreachable. Consequently, the classic phase requires a minimal Kg^0.027, leading to the 
observed reentrant phenomenon. 

The morphology of the crumpled wires can be quantified by the number of loops, the 
distribution of loop sizes and contact points, or the fractal dimension of the structure. In 
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FIG. 2: Morphological phase diagram of the spiral/plastic/classic morphologies obtained numeri- 
cally for /=50, where the control parameters are the static wire friction fj, s t and the elastic yield 
curvature kq. Each point is averaged over 9 realizations. 



this letter we focus on the scaling of the number of loops Ni as a function of the dimensionless 
packing density <fi = d- L/ (vrr 2 ), where L is the inserted wire length. A loop is denned as an 
area that is surrounded by a wire segment with only one inner contact of inward wire surfaces. 
In experiments, the number of loops was counted and the total wire length was obtained by 
standard image analysis of digital images that were taken during the experiment. Gomes 
et al. 8| report a pronounced shoulder for 0<O.O32, a power law asymptotic dependence 
Ni ~ cp with 7=1.8±0.2 and a maximum packing density of 0=O.14±O.OO6. The results 
of our simulations and experiments are given in Figj3l We find exponent 7 to vary slightly 
for different morphologies. For the classical morphology we measure 7=1.75±0.03, for the 
plastic morphology 7=1.85±0.04 and for the spiral 7=0 up to 0=1, as only two loops are 
present. Also a universal maximum packing density could not be observed. 

For a macroscopic analysis of the packing process, we measured the stiffness of the crum- 
pled structures as function of the packing density <fr via the force acting on the nodes in the 
insertion channels. As before, the simulation is stopped when the system is jammed, i.e. 
when straight parallel force lines from contacting wire segments exist between the injection 
channels (see FigHJ). 

For computational reasons, these measurements were performed for different system sizes 
up to /=35. Parameters were chosen as kq=0.08, s=0.06 and n st =5, corresponding to the 
classical phase. We find the force to scale as a power-law of the form 

F oc (cf> c - 4>f (2) 




FIG. 3: Scaling of the number of loops for the classic (/=50) and plastic (/=100) morphologies. To 
compare different systems sizes, L/Lq was used as packing ratio (total inserted wire length L and 
initial one Lq=2R). The number of loops for both morphologies follows the power law iVj ~ (/)"'. 




FIG. 4: The jammed classic morphology. Force chains formed by contacts are clearly visible. Con- 
tact forces are represented by the thickness of the line segments, while the grey scale of the wire 
backbone represents its local elastic bending energy. Since this system includes plastic deforma- 
tions, highest curvature does not necessarily imply highest bending energy. 



with exponent /3=-1.43±0.02. In FigfSl </> c =0.46±0.01 which is substantially smaller than 
the theoretical limit C =1. For the purely elastic case of the classic morphology, we find 
/3=-2.05±0.02 and C =O.54±O.O2 (data not shown). Note that simulations for the spiral 
phase exhibit <p c ~ 1 with a deviant force scaling relation, in excellent agreement with the 
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FIG. 5: The insertion force F as function of 4>. Simulation parameters are &i=1000, /c2=60, /=25, 
Ke=0.08, /i s t=5, s=0.05. The best fit for the divergence is a power law F oc {(j) c — (j)f with /3=-1.43 
(solid line). 

analytical results jsj]. 

In this letter we presented three morphologies into which wires crumple inside a circular 
cavity. Using a discrete element model that incorporates plastic deformation and static 
friction, we showed that plasticity and friction are the two essential parameters determining 
the morphology. Three different experimentally found morphologies could be reproduced 
in silico by matching wire materials and friction, and the associated phase diagram was 
constructed. We showed the existence of a reentrant phenomenon due to two mechanisms 
for the stabilization of structures, friction and plastic deformation. The crumpled structures 
were analyzed by counting the number of loops as function of the packing density, showing 
a power-law behavior for experiments and simulations. The insertion force, which is difficult 
to obtain in experiments, was determined numerically and exhibits a power-law divergence 
with critical packing fractions substantially lower than 1. The same power-law divergence 
(with C ^O.75, /3~-1.85 [h], [jjjl) was also found for the forced crumpling of membranes in 



151 ] . In wires, 



three dimensions, although the divergence in membranes comes from kinks 
on the other hand we have no kinks and it is most probably a consequence of the decreasing 
size of the newly generated loops. 

The presented work demonstrated the crucial role of plasticity and friction in crumpling 
processes, with its characteristic morphology. It is natural to generalize this work to the 



forced crumpling of membranes or the packing of DNA into viral capsids 



16 



12]. Macro- 
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scopically, the role of the system size is not fully understood and certainly needs clarification, 
also with regard to the morphological phase diagram. Furthermore, it is still an open ques- 
tion how other geometrical parameters such as the injection angle and cavity shape influence 
the morphological phase diagram. 

This work was supported by the Grant TH-06 07-3 of the Swiss Federal Institute of 
Technology Zurich. 
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